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ABSTRACT 

Model atmospheres have been computed for M dwarfs that are strongly irradiated by nearby hot 
companions. A variety of primary and secondary spectral types are explored in addition to models 
specific to four known systems: GD 245, NN Ser, AA Dor, and UU Sge. This work demonstrates 
that a dramatic temperature inversion is possible on at least one hemisphere of an irradiated M 
dwarf and the emergent spectrum will be significantly different from an isolated M dwarf or a black 
body flux distribution. For the first time, synthetic spectra suitable for direct comparison to high- 
resolution observations of irradiated M dwarfs in non-mass transferring post-common envelope binaries 
are presented. The effects of departures from local thermodynamic equilibrium on the Balmer line 
profiles are also discussed. 
Subject headings: stars: atmospheres, binaries: close, radiative transfer 



1. INTRODUCTION 

Synthetic spectra from model atmospheres are fre- 
quently used in the analysis of observed spectroscopic 
and photometric data. For the most part, the models 
are sufficiently detailed to test the current theoretical 
understanding of stellar and sub-stellar mass objects at 
various stages in their evolution. However, the major- 
ity of model atmospheres are intended for comparisons 
to isolated stars and are not appropriate for many short 
period binaries. A number of binary systems have or- 
bital separations small enough so that one of the binary 
members is significantly heated by its companion. In or- 
der for synthetic spectra to be useful in such cases, the 
standard "isolated" modeling approach must be replaced 
by one that includes the effects of irradiation. 

Post-common envelope binaries (PCEBs) are examples 
of systems where large effects due to irradiation have 
been observed (e.g. Z Cha, Wade & Home (1988); GD 
245, Schmidt et al. (1995); GD 444, Marsh & Duck 
(1996)). Many PCEBs are cataclysmic variables (CVs) 
with secondaries that, in addition to being irradiated, 
have over filled their Roche lobe and transfer mass to the 
primary via accretion streams or a disk. CVs, while very 
interesting, are not necessarily the best choice for study- 
ing the effects of irradiation on cool stars. PCEBs that 
do not have on-going mass exchange offer several advan- 
tages. The absence of accretion ensures that the primary 
is the only major source of external heating for the sec- 
ondary. CVs have additional sources of radiation energy 
(e.g., from the accretion disk) that are complicated to 
model and could lead to shadows or bright spots on the 
secondary. An additional advantage of PCEBs without 



mass exchange is that the secondaries have not over filled 
their Roche lobe and are therefore more spherical. Most 
non-mass transferring PCEBs are labeled pre-CVs and 
often contain a main sequence (MS) star in close orbit 
around a much hotter white dwarf (WD) or sub-dwarf 
(sdOB). The current study will be restricted to pre-CVs 
with orbital periods less than 16 days or orbital sepa- 
rations less than about 3Rq (i.e., those t hat could pos- 
sibly become a CV within a Hubble time; iHillwig et alJ 
(20001)). For an excellent re view of detac hed binaries 
containing WD primaries, see l)Marsh "2000") . 

The average effective temperature of the primary (the 
WD or sdOB) in pre-CVs is about 50,000K with a few 
reaching 100,000K. At just a few R© away, the much 
cooler secondary star is significantly heated by the pri- 
mary's radiation. Therefore, the secondary's irradiated 
atmosphere is regulated by both extrinsic radiation from 
the primary and intrinsic energy supplied by internal 
nuclear reactions (or, in the case of a brown dwarf sec- 
ondary, leftover gravitational energy from formation). 
Depending on the temperatures, pressures, and chemical 
composition of the secondary's atmosphere, a fraction 
of the extrinsic radiation will be reflected via scattering 
by grains, molecules or electrons (or a combination of 
these) while the remaining is absorbed and re-radiated. 
Furthermore, the two sources of energy may have spec- 
tral energy distributions (SEDs) that peak at very dif- 
ferent wavelengths which implies that a broader category 
of opacity sources (spanning the EUV to the far IR) be- 
come important in shaping the seconday's atmospheric 
structure and emergent spectrum. The situation is also 
complicated by the fact that only one hemisphere of the 
irradiated companion is heated at any given time. 
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Fig. 1. — An illustration of the irradiation model. The primary 
(left) is treated as a point source at distance a from the center 
of the secondary (right). The secondary is divided into two hemi- 
spheres: irradiated and non-irradiated. The atmosphere on either 
hemisphere is modeled as a collection of concentric shells with the 
height of the atmosphere equal to R2 - Rl. 



Only a few theoretical studies applicable to the at- 
mospheres and spectra of pre-CVs have been published. 
One of the more recent works investigated the effects 
of irradiation in th e eclipsing binary BE U rsae Ma- 
joris using CLOUDY IjFerguson fc James! I1994D . A nar- 
row range of M dwarfs (3400K < T^s < 3800K) lo- 
cated near a 10,000 K bla ck body have also been mod- 
eled (Brett & Smithl ll99.1 here after BS93). A few ear- 
lier papers provided critical insight into the problem of 
irradiatio n but did not consider situations relevant fo r 
pre-CVs HVaz fc Nordhmdiri985tlNordhiTid fc VaH[T99?1). 
In addition, the necessary opacity data and spectral line 
databases needed to compute realistic synthetic spec- 
tra were not available at the time. Consequently, no 
synthetic spectra have been published that are detailed 
enough to provide useful comparisons to high-resolution, 
phase-resolved, observations of pre-CVs. Furthermore, 
most light curve models continue to use black body SEDs 
or, at best, non-irradiated atmosphere models. 

In the following sections, models for a variety of ir- 
radiated pre-CV secondaries will be presented. The 
changes in the atmospheric structure, the chemical com- 
position, and the spectra will be described. Detailed, 
line-blanketed, synthetic spectra will also be presented 
for a few specific pre-CV systems. 

2. MODEL CONSTRUCTION 

All calculations presented here were pro- 
duced using the PHOENIX model atmosphere code 
llHauschildt & Baron..l9991 PHOENIX includes an equa- 
tion of state and radiative transfer solver th at is suitable 
for modeling a broad range of both hot IjAufdenberd 
HqOI) and cool (AUard ct al. 2001) objects across the 
H-R diagram. This flexibihty makes PHOENIX well 
suited for modeling systems that contain two extremely 
different objects (e.g. M dwarfs and WDs). 

The secondary may be conceptually divided into two 
distinct hemispheres; irradiated and non-irradiated (see 
Figure^. The irradiated hemisphere receives incoming 
radiation over all angles between 0° and 90° with re- 
spect to the surface normal. The effects of heating by 
the primary will vary across the surface with decreased 
heating near the boundary between the two hemispheres. 
In order to simplify the simulations, each hemisphere was 
modeled separately using a 1-D, spherically symmetric, 
atmosphere designed to represent the average thermal 



and spectroscopic properties. Since pre-CVs are non- 
mass transferring systems, the secondary does not over- 
flow its Roche lobe and, hence, should be reasonably well 
approximated by a sphere. Consequently, the same mass 
and effective gravity (set at a reference optical depth) was 
used for both hemisphere models. However, despite hav- 
ing the same bulk properties, maintaining a consistent 
temperature structure and entropy {S) at the bottom 
of both atmosphere models for the irradiated and non- 
irradiated hemispheres is not guaranteed. Differences be- 
tween the temperature structures at the bottom of each 
model can exist because the effects of horizontal energy 
flow have not been included and the two hemispheres 
have been essentially decoupled. While horizontal en- 
ergy flow is very likely to be present, including this in the 
simulations would require the solution of the full, multi- 
D, radiation-hydrodynamical problem that is currently 
beyond the scope of this paper. The standard approach 
for dealing with horizontal energy flow is to adjust the 
effective temperature of one model so the entropy at the 
deepest layer is the same for both hemispheres (Vaz & 
Nordlund 1985; BS93). The details of this problem and 
how it is handled are discussed in a subsequent section. 

All models were computed assuming hydrostatic and 
radiative-convective equilibrium. Convection was in- 
cluded using the standard mixing length theory with 
the mixing length parameter set equal to two. The at- 
mospheres were divided into 100 concentric shells and 
the spherically symmetric radiative transfer equation was 
solved while explicitly including the incident, frequency 
dependent, radiation field from the primary. Spherical 
symmetry, as opposed to the more traditional plane par- 
allel geometry, was chosen because it provides a more 
realistic description of the incoming radiation at graz- 
ing angles that should not penetrate too deeply into 
the atmosphere. At some angles, the radiation can pass 
completely through the top of the atmosphere. The de- 
tails of the radiative transfer solution may be found in 
iHauschildt fc BaronI l)1999|) (and r efs. therein) a r id de- 
tails of the irradiation are found in 'B arman et alJ l)200lL 
,2002); however, for clarity, a few of the important points 
are repeated. 

A secondary with radius Rg intercepts an incident lu- 
minosity from a primary of radius Rp at distance d with 

flux Fp given by ttR-^ ( if ) -^p- This absorbed energy 

must be distributed over and re-radiated by some frac- 
tion of the full secondary surface area {AnRf). The 
three most commonly explored physical scenarios are re- 
radiation by the full surface area, by the irradiated hemi- 
sphere only, or by a single point. Since the models must 
conserve energy, the incoming luminosity must balance 
the out-going luminosity of re-radiated extrinsic energy. 

In which case, ttR^ i-j-) Fp — Area x F^ 

further assumes that Fmc 
incident fluxes can be expressed as, 

Fu.cx^af^'] Fp^^ 



out- If one 
-Fout, the monochromatic 



(1) 



where Fp\ are the monochromatic fluxes from the pri- 
mary surface, and d is the distance from the primary 



surface to the secondary surface (note, d 



Rp — Rs, 



where a is the orbital separation) . The parameter a sets 
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Fig. 2. — Temperature versus pressure for irradiated (top) and 
non-irradiated (bottom) atmosphere models. Solid lines are the 
results of PHOENIX simulations while dashed lines are taken from 
Fig. 3 of Brett & Smith (1993). Filled circles indicate the location 
of the radiativc-convective boundary in the PHOENIX models. The 
differences between the two sets of calculation arc most likely due 
to the different treatments of the equation of state and opacities. 



the redistribution of the absorbed incident energy over 
the secondary's surface. All incident energy being ab- 
sorbed and re-radiated by the heated hemisphere corre- 
sponds to a = 0.5. Re-radiation by the entire secondary 
surface corresponds to a = 0.25 and a = 1.0 is for re- 
radiation by a single point on the surface. For a more 
detailed description of the energy balance in an irradi- 
ated binary compan ion and th e deve lopment of a simi- 
lar a parameter, see iPaczvnskil l|1980() . Unless otherwise 
stated, the models presented below were calculated with 
OL = 0.5. This choice is motivated by the fact that many 
short period binaries have tidally locked secondaries. 

In all cases, the incident radiation was approximated 
by an isotropic radiation field consisting of 64 in-coming 
intensities (corresponding to 64 incident angles) at the 
top most layer (ti.2 = 0, where ri.2 refers to a reference 
optical depth measured at 1.2/im) of the model. Isotropic 
illumination (instead of heating at a fixed angle) was 
chosen because it better approximates the situation il- 
lustrated in Fig. n While this is clearly a simplification 
of the true incident radiation, isotropy roughly accounts 
for the fact that the heated hemisphere receives radiation 
from the primary at all angles between and 90 degrees 
(with respect to the surface normal). Even when the inci- 
dent radiation is assumed to be isotropic, there are still 
benefits to using spherical geometry. Spherical geome- 
try, combined with isotropic irradiation, accounts for the 
preferential heating of the upper atmospheric layers com- 
pared to the deeper layers for regions near the boundary 
between the two hemispheres. Across the heated hemi- 
sphere, the upper layers always receive extrinsic radia- 
tion while the deeper layers experience considerably less 
heating near the boundary compared to the substellar 
point. By virtue of the translational symmetry inherent 
in plane parallel models, all inward directed intensities 
are propagating toward the center of the object. Thus, 
for some situations, plane parallel geometry could overes- 
timate the heating of deeper atmospheric layers and un- 
derestimate the heating of upper layers. The differences 
between plane parallel and spherical models are maxi- 



mized when the atmospheric extension is large compared 
to the mean free path of a photon (e.g., in low gravity 
objects). For a spherical model, a large extension has the 
effect of decreasing the range of angles over which inci- 
dent intensities actually propagate toward the bottom of 
the photosphere. In the future, the heated hemisphere 
will be modeled as a collection of small annular regions 
where the extrinsic radiation is incident along a single an- 
gle (or a narrow range). In this case, the advantages of 
spherical geometry over plane parallel will be even more 
significant. 

Since pre-CV primaries have spectra that are very dif- 
ferent from black body SEDs, the incident intensities 
were taken from synthetic WD spectra (also calculated 
with PHOENIX; Barman et al. 2000). In this way, a more 
realistic, wavelength dependent spectrum was used for 
the incoming radiation. Unless otherwise stated, the WD 
models have \og{g) — 8.0 (cgs units), 1 M©, and 10~^ so- 
lar metal abundances. When modeling a secondary for a 
specific pre-CV, a new model atmosphere and synthetic 
spectrum were used that closely match the primary spec- 
tral type and, if available, the observed fluxes. 

The model calculations for both the primary and the 
secondary included more than 50,000 wavelength points 
between lOA and 1000/im with a resolution of about lA 
from the UV to the near-IR (this resolution may be in- 
creased as necessary). The equation of state and opacity 
setup was similar to that us ed in the Next Gen model at- 
mosphere grid of Hauschild t et al.l (|1999D but has since 
been updated to include more recent chemical and opac- 
ity data ( AUard et al. 2001). Local thermodynamic equi- 
librium (LTE) was assumed for the majority of the calcu- 
lations; however, the effects of non-LTE are explored for 
the hydrogen atom. All models for the secondaries used 
a solar composition that included 40 of the most impor- 
tant atomic elements from H to La (atomic numbers 1 
through 57) as well as their important ions. The atomic 
data for the ener gy levels and bound-bound transitions 
are from iKuruc J 



J1994I) and lKurucz &BeIil llT995l). The 

molecular opacities include HtO IjPartridee k, Schwenkd 
1997), TiO (Schwcnkc 1998), VO (hnelist provided by R. 
Free dman, 2001, p rivate communication) all diatomics 
from lKuruc3 l|1993j), and all lines from the HITRAN and 
GEIS A databases l|Rothman et all Il992t IHusson et alJ 
I1992J) . Colhsion induced absorption (CIA) opacities are 
also include d for Ha, N2, A r, CH4, and CO2 according to 
(jGruszka fc Borvsowlll997l and references therein). The 
total number of atomic and molecular lines currently 
available in PHOENIX (version 13) is ~ 700 milhon. Scat- 
tering (Thomson, Rayleigh, and Mie) was also included 
and assumed to be isotropic. 

The PHOENIX model atmosphere code uses the effec- 
tive temperature as an input parameter that specifies 
(through the Stefan-Boltzmann formula) the net flux (or 
luminosity) that the model will have. To avoid confusion 
between the effective temperature of irradiated and non- 
irradiated models, the variable Tog' will be used only to 
refer to the effective temperature of non-irradiated mod- 
els. In the case of irradiated models, Tint will be used 
to describe the intrinsic effective temperature; i.e. the 
effective temperature the model would have if irradiation 
was not present. 

3. RESULTS 
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Fig. 3. — Temperature as a function of optical depth at 1.2/im 
for irradiated (solid lines) and non-irradiated (dashed lines) at- 
mospheres. The irradiated atmospheres have Tp = 20,000K to 
lOO.OOOK (from top to bottom) in steps of 10,000K. For all irra- 
diated models, the orbital separation is 1 R0 and Tint = 3000K. 
The non-irradiated models have T<,ff = 5600K, 4000K, and 3000K 
(from top to bottom). Filled circles indicate the location of the 
radiative-convective boundary. 



3.1. Comparison to the BS93 Results 

Figure [21 compares the temperature-pressure (T-P) 
profiles of non-irradiated and irradiated models to simi- 
lar calculations from BS93 (see their Fig. 3). The mod- 
els have Tint = 3400K, log(g) = 5.0 (cgs units) and Ms 
= 0.22 M0 (note, BS93 used plane parallel models and 
so Ms is irrelevant). The PHOENIX irradiated model in- 
cluded a primary located at an orbital separation of 0.39 
Rq with Tp = 10,000 K, Rp= O.O319R0, and a = 1.0. 
These parameters were chosen to match the incident flux 
used by BS93. 

The models in Fig. |21 illustrate the basic properties 
of an irradiated atmosphere. As should be expected, 
substantial heating occurs in the upper layers of the ir- 
radiated atmosphere (low -Pgas) where temperatures ex- 
ceed those in the non-irradiated case by 3000K or more. 
The temperature also increases throughout the atmo- 
sphere, even at regions with high gas pressures, result- 
ing in a different inner adiabat compared to the non- 
irradiated model (the reasons for this are discussed be- 
low). Also, the radiative-convective boundary nearly 
reaches Pgas ~ 10^ dynes cm^^ as convection retreats 
to deeper layers in the irradiated atmosphere. 

Both PHOENIX and BS93 models predict a similar tem- 
perature inversion at the top of the irradiated atmo- 
sphere and similar adiabats at high Pgas- However, the 
BS93 profiles are noticeably cooler than the PHOENIX pro- 
files at nearly all Pgas . Above the temperature minimum 
in the irradiated case (near Pgas — 10^ dynes cm^^), the 
temperatures are more than 800K higher in the PHOENIX 
models. These differences are likely due to the use of 
"straight mean" (SM) molecular opacities and a lack 
of any atomic line opacity in the BS93 models. The 
lack of atomic line opacity can explain the differences 
at high pressures due to an underestimated backwarm- 
ing affect. At low pressures, it is difficult to pinpoint 
the exact causes of the differences but, as pointed out 
by BS93, SM opacities tended to overblanket their mod- 
els which could have caused preferential heating of the 



outer layers. Note that, when the BS93 models were 
produced, reliable atomic and molecular line-lists were 
only just becoming available and thus the use of SM 
opacities was quite common. Large line databases, how- 
ever, have been steadily produced over the past few years 
and are incorporated into PHOENIX using a direct opacity 
sampling method. A second major difference between 
the two sets of models is the completeness of the equa- 
tion of state. PHOENIX includes a far greater number of 
species (literally many hundreds of molecules) in the so- 
lution of the chemical equilibrium equations compared 
to the BS93 models. Differences in the chemical equi- 
librium can lead to very different concentrations of pri- 
mary opacity sources and ultimately very different tem- 
perature structures. Despite the differences between the 
models, the fact that two independent atmosphere codes 
produce qualitatively consistent results is encouraging. 

3.2. Atmospheric Structures 

To explore the effects of irradiation by primaries with 
different Tp values, a sequence of models was calculated 
for a secondary with Tint = 3000K, log(g) = 5.0, Ms= 
0.22 M0 and solar metalicities. The primary was a 
O.OI4R0, log(5) = 8.0 (sub-solar metahcity) WD located 
IR0 from the secondary. Tp was increased from 20,000K 
to 100,000K by steps of 10,000K. With this setup, the 
incident flux spans almost 3 orders of magnitude with 

-0.7 < l0g(Pinc/Pi„t) < 2.1. 

The temperatures as functions of optical depth (T1.2) 
are shown in Fig. for each irradiated model atmo- 
sphere. The models are fully converged (i.e., energy is 
conserved to better than 2% at all layers) and the atmo- 
spheric chemistry and opacities are consistent with the 
changes in temperature. For the lowest primary tempera- 
ture (Tp = 20, OOOK), a significant temperature inversion 
forms in the upper regions of the atmosphere while the 
deeper layers (near T1.2 = 1) remain fairly similar to a 
non- irradiated atmosphere with Teff = Tint- As Tp in- 
creases, the temperature inversion also increases and the 
temperature minimum moves to deeper layers. With Tp 
= 100, OOOK, the temperatures at the top of the atmo- 
sphere reach 19, OOOK and the T-minimum is located at 
T1.2 < 1.0. Once Tp is higher than 30, OOOK, the inner 
temperature profile is no longer similar to that of a stan- 
dard M dwarf, despite the fact that Tint was held fixed 
for this sequence of models. 

The large range of atmospheric temperatures and pres- 
sures displayed in Fig. |31 is accompanied by a very di- 
verse chemical make-up. In non-irradiated M dwarfs, 
molecules (H2, H2O, CO, TiO, etc.) are the major 
sources of opacity throughout the atmosphere. For an 
irradiated atmosphere with Tp = 20, OOOK, the incident 
flux raises the temperatures in the upper atmosphere so 
much that most molecules cannot survive. Figure 0] illus- 
trates the changes to the number densities of H2, H, H+, 
He, and He+. When Tp = 20, OOOK, the total concentra- 
tion of molecular species (Nmois/Ngas) at ri.2 < 10~^ de- 
creases by more than 5 orders of magnitude compared to 
the non-irradiated case. Also, a significant fraction of hy- 
drogen is ionized near the top of the atmosphere. When 
Tp > 40, OOOK, molecules are no longer important in the 
photosphere and the concentration of H2 at the temper- 
ature minimum is less than 10^^. Other molecules (e.g. 
CH4 and CO) follow a similar pattern as H2 but with 
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Fig. 4. — The concentrations (number density A'^ divided by the total number density Ntot) of H, H+, H2, He, and He+ are shown 
for a Teff = 3000K non-irradiated model and the irradiated models with Tp = 20,000K, 40,000K, 60,000K, 80,000K, and 100,000K. The 
corresponding temperature structures are shown in the bottom right panel. 



smaller concentrations at all layers. At even larger Tp, 
the atmosphere becomes mostly composed of atoms and 
ions and, when Tp = 100, OOOK, hydrogen is almost com- 
pletely ionized down to T1.2 = 0.1 and He is ionized down 
to T1.2 — 10^^. The ionization in the upper half of the 
atmosphere dramatically increases the number density of 
free electrons which, in turn, causes the number density 
of H+ and He+ to increase. The heating of the inner 
atmospheric layers also leads to increased electron den- 
sities and concentrations of H+ and He+. Overall, the 
number densities of electrons, H+, and He+ follow the 
same pattern as the temperature structure. Note that, 
while the solution of the EOS in PHOENIX allows for the 
possible formation of dust grains and other condensates 
(as described in AUard et al. 2001), all of the models 
presented here have atmospheric temperatures that are 
too hot for grains to form. For intrinsically cooler sec- 



ondaries (e.g., brown dwarfs and extrasolar planets) and 
relatively low irradiation, grain formation will be very 
important (,Bar man et al. 2001). 

3.3. Entropy Matching 

The increase in temperature at large depths (i.e., large 
Pgas) and throughout the convection zone has been dis- 
cussed in previous papers (Vaz & Nordlund 1985; BS93). 
As described in section |21 the model atmospheres are 
constrained to be in thermal equilibrium and, therefore, 
each layer in the convection zone was required to trans- 
port a constant total energy. Furthermore, in a purely 
adiabatic convection zone, the energy transport is deter- 
mined by the adiabatic temperature gradient. Thus, a 
temperature increase near the top of the convection zone 
must lead to an increase at all other layers in the convec- 
tion zone to ensure a constant luminosity is maintained 
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Fig. 5. — The bolometric albedo for the irradiated atmospheres 
in Fig. |21(filled symbols). Error bars that assume a 2% error in the 
net flux of the model are given for each Tint = 3000K model. For 
large i^inci the error bar s are smaller than the s ymbol size. Also 
plotted are albedos from INordlund &: Vad 119901) (diamonds) and 
BS93 (asterisk) along with three PHOENIX models (pluses) specifi- 
cally for comparing with these earlier works. 



at all layers of the atmosphere. The now well known con- 
sequence of these constraints is that convective irradiated 
and non-irradiated models with the same effective tem- 
perature (i.e. Tint = Tes) do not reach the same adiabat 
at depth (as illustrated in Figs. 13 and 13). Furthermore, 
since the material below the photosphere is likely to be 
well mixed, the irradiated and non-irradiated model at- 
mospheres should have the same entropy at depth if they 
are to represent the two different hemispheres of the same 
object (Vaz & Nordlund 1985 ; BS93). In most cases, en- 
tropy matching can be achieved by lowering Ti„t of the 
irradiated model until the T-P profile matches the non- 
irradiated model in the convection zone. The difference 
between Tint and Tcs for the matching irradiated and 
non-irradiated atmospheres is a measure of the bolomet- 
ric reflection albedo and the energy that is horizontally 
redistributed in the photosphere (BS93). 

The entropy matching technique is designed to mimic 
the anticipated conditions at large optical depths within 
the convection zone. These layers should undergo effi- 
cient mixing over very short timescales ensuring a con- 
stant internal entropy across both hemispheres. Also, 
the decrease in the atmospheric temperature gradient on 
the heated hemisphere will slow the rate at which inter- 
nal energy can escape from that side of the secondary 
implying a smaller intrinsic luminosity on the heated 
hemisphere compared to the non-heated side. This is 
consistent with the fact that entropy matching models 
will always have intrinsic effective temperatures that are 
less than (or equal to) the Tcs on the non-irradiated side. 
Entropy matching, however, cannot deal with the poten- 
tial horizontal flow of energy in the radiative zones. 

To flnd pairs of entropy matching models, a sequence 
of non-irradiated models was computed with 2800K < 
Toff < 6000K. Then, using the entropy at the bottom of 
an irradiated model, the matching effective temperature 
(Tinatchcd) was found by linearly interpolating along a 
S'(Toff ) curve at the same Pgas in the non-irradiated mod- 
els. ThisjnethodjWhichJs very similar to that described 
bv IVaz fc Nordlundl l)1985j) . was applied to each of the 



models shown in Fig. |21 The matched, non-irradiated, 
atmospheric structures for the Tp = 20,000K, 50,000K, 
and 100,000K models are also shown as dashed lines and 
have, respectively. Toff = 3041K,3857K, and 5910K. Re- 
call that the corresponding irradiated atmosphere models 
all have Tjnt = 3000K. 

The two effective temperature values for the pair of en- 
tropy matching models can be used to obtain a value for 
the bolometric reflection albedo (w) via the expression. 



(1 - w)Fi, 



^ (^matched ^off) 



(2) 



(BS93). The albedos versus the ratio of incident to in- 
trinsic fluxes are shown in Fig. For very large incident 
flux the albedo approaches 1. This is expected since the 
atmospheres become radiative down to large depths as 
Tine increases. From Fig. the radiative convective 
boundary is at T1.2 ^ 1 when logTino/Tjnt > 0.5. Atmo- 
spheres that are completely in radiative eq uilibrium have 
an albedo exactly equal to 1.0 l)Eddingtorl l926'). For low 
incident flux the albedo decreases but is still well above 
the canonical value of 0.5 for convective atmospheres. 
In general, the PHOENIX albedos are consistent with the 
recently published albedos for a v ariety of mod els with 
Tint > 4000K and logTi„o/Ti„t < 1 (|C]aTeH2nmL Fig. 4). 

All models in this paper were considered converged 
when the flux was conserved to 2% or better at all lay- 
ers. At this level of convergence, no appreciable improve- 
ments in the models are obtained by further iterations. 
None the less, a small error exists in the net flux of 
the models which translates into a small error associated 
with T^^i:,tchcd ill ^l- El For small values of Tine, the dif- 
ference between Tmatched and Tefj is also small (^ 50K) 
which means that even a 2% error in the flux can result 
in a large error in w. Error bars corresponding to a 2% 
error in the model's net flux are shown for each albedo 
(Tint — 3000K models) in Fig. |S| Since the error in w is 
also proportional to 1/Tinc, the error bars decrease with 
increasing Tine (and the difference between Tmatehed and 
Toff is much larger). 

Comparison PHOENIX models were also constructed 
to match the inp ut parameters of several models from 
INordlund fc Vad <)1990(1 (see their Table 2; non-grey 
models with mixing length parameter equal to 2) in ad- 
dition to the BS93 models shown in Fig. The PHOENIX 
abledos for these models arc plotted in Fig. El and agree 
wefl with Nordlund & Vaz (1990) while the BS93 albedo 
is less than half the value for any of the others. It is 
difficult to determine the exact reasons why the BS93 
result is so much lower than the PHOENIX albedo. In or- 
der to obtain an albedo of 0.2 for the PHOENIX model, 
Tniatehod would nccd to bc '^ 3800K, instead of the cur- 
rent value of -- 3600K. Note that the BS93 models 
suffered from many competing simplifications that in- 
cluded SM opacities, the absence of atomic fine opac- 
ity, and the use of underestimated astrophysical values 
for the TiO oscillator strengths. The impact of these 
simplifications can be seen by comparing the BS93 non- 
irradiated models to Brett's later, much improved, non- 
irradiated M dwarf models fBrett 199^. The ten ipera- 
tureat the bottom of a Toff = 3800K model from iBrettI 
|l995| is nearly lOOOK hotter than the bottom of the 
non-irradiated model shown in BS93 with the same Toff 
(their Fig. 7). However, at Toff = 3400K, the differences 
at the bottom are on the order of a few lOOK. Using the 
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Fig. 6.— Spectra for the irradiated models shown in Fig. 01 Tp = 20,000K to 50,000K in the lower panel and Tp = 60,000K to 100,000K 
in the upper panel (with Tp increase from top to bottom). To facilitate comparison, the resolution was reduced to 15A and alternating line 
styles were used. The dashed line is a best fitting black body spectrum (with T^s = 4840/C) for the Tp = 50,000K model.. 



improved iBretd IjlQQIif l non-irradiated models to find the 
Tlnatchod for the BS93 irradiated model leads to an albedo 
very comparable to the PHOENIX value. The good agree- 
ment in Fig. Elbetween the PHOENIX and BS93 irradiated 
models suggests that, for this particular set of parame- 
ters, the extrinsic heating offsets many of the opacity 
limitations in the BS93 models. 

3.4. Synthetic Spectra 

The changes to the temperature structure and chemi- 
cal composition have dramatic consequences for the spec- 
tra emerging from the illuminated hemisphere. The syn- 
thetic spectra for each of the structures in Fig. O are 
displayed in Fig. |^ The atmospheric model irradiated 
by the coolest WD (Tp = 20, OOOK) retains many of 
the infrared spectral features common to M dwarfs, e.g. 
TiO and VO absorption bands. However, a substantial 



amount of flux emerges at short wavelengths with a dense 
spectrum of narrow emission lines. As Tp increases, the 
flux emerges from two fairly distinct spectral regions di- 
vided by a pronounced Balmer-edge at 3647A. When Tp 
reaches 30, OOOK, the broad spectral features in the red 
portion of the secondary's spectrum disappear and are 
replaced by a few narrow emission lines. The blue-UV 
portions of the spectra, however, have complex mixtures 
of emission and absorption lines. The strongest of these 
hues are from Fe I, Fe II, Mg I, Mg II, and Ni I. When 
Tp — 100, OOOK, many additional emission lines from 
Fe III, Al II, Ti II, Ni II, and Si I appear in the UV- 
blue spectrum. 

This division of the spectrum into two distinct com- 
ponents is a consequence of the two very different tem- 
perature regimes in the atmosphere. The hot, inverted, 
outer temperature structure produces the blue-UV por- 
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Fig. 7. — Right panel: the ratio of departure coefficients for the levels involving Hq, H^, and H-y as functions of optical depth at 1.2/im. 
Each ratio is proportional to the non-LTE line source function. A value other than 1 indicates departure from LTE. Left panel: comparison 
between LTE (dashed) and non-LTE (solid) line profiles for the corresponding Balmer lines. Each panel is 20A wide. 



tion of the spectrum. Note, r = 1 at 3000A is at roughly 
Ti.2 — lO"**. The red part of the spectrum forms near 
Ti.2 = 1, while the core of most emission lines form at 
smaller ti.2. For all of the cases presented in Fig. the 
incident radiation is sufficient to produce Balmer, He, 
and metal emission lines. 

A best fitting black body SED, with Tcff = 4840K, is 
also shown in Fig. Elfor A > 4500A and Tp = 50,000K. 
Apart from the emission lines, the red portions of the 
spectra resemble black body SEDs when Tp > 30,000K. 
This is mostly due to the roughly isothermal T — P struc- 
ture near the temperature minimum that occurs close to 
continuum forming region (ti.2 ~ 1.0). Spectral features 
diminish where the temperature structure is flat because 
the LTE source function cannot vary with depth. At 
shorter wavelengths, the best fitting black body underes- 
timates the emergent flux by orders of magnitude. Also, 
for Tp > 60, COOK, a steep Blamer-edge forms in emis- 
sion. A steep Blamer-edge could explain why, for certain 
binaries, it is difficult to obtain a single effective temper- 
ature for the secondary when using black body SEDs to 
fit ubv photometry. 

3.5. Hydrogen in non-LTE 

A striking feature of the spectra shown in Fig. are 
the hydrogen emission lines. However, one caveat should 
be mentioned for the model results discussed so far; in 
all of the above cases LTE was assumed. Given the non- 
local origin of a primary source of energy (the incident 
flux), the conditions in the upper atmosphere, where gas 
pressures are low and temperatures are high, are well 
suited for significant departures form LTE. To test the 
importance of non-LTE effects, the Tp = 100,000K model 
was recalculated with neutral Hydrogen allowed to de- 
part from LTE. For this particular test, the Hydrogen 
model atom included 31 levels, and 435 primary transi- 
tions (all bound-bound transitions with \og{gf) > —2.0) 
were included in the solution of the statistical equilib- 
rium equations. For details of th e numerical methods 
used to solve the rate equations, see lHatischildt fc BaronI 

Earn. 



The LTE and non-LTE line profiles for Hq, Rp, and 



H.y are compared in Figure [7| Both LTE and non-LTE 
cases have Balmer wings in emission while dramatic dif- 
ferences exist between the line cores. The non-LTE line 
profiles all have reversed cores while the LTE lines have 
cores that are substantially brighter than the non-LTE 
case. The non-LTE Hq, line core drops slightly below the 
continuum flux level while the cores of H^, and H-y are 
only slightly above the continuum. The ratio the LTE 
and non-LTE Balmer line fluxes, for this case, are not 
significantly different. The non-LTE flux ratio Hq,/H^ is 
^ 0.52 compared to ~ 0.57 for the LTE case. The non- 
LTE flux ratios for Hq/H^ and H^/H.^ are ^ 0.41 and 
~ 0.80, respectively. 

In both non-LTE and LTE cases, the line wings form 
near ri.2 = 1 while the core of the lines form much higher 
in the atmosphere near T1.2 = 10^"*. In LTE, the depth 
dependent line source function is simply a black body 
with temperature equal to the local gas temperature. 
Consequently, the LTE line profile mimics the temper- 
ature inversion (top curve in Fig. O which changes by 
~ 9000K between these two optical depths. However, 
when the LTE approximation is removed, noticeable de- 
partures from the LTE atomic level populations occur for 
the n = 2 level, which is over populated by a factor 10 
(compared to the LTE populations) at T1.2 = 10""*. The 
n = I and 3 < n < 6 levels are also slightly over popu- 
lated in the upper atmosphere, but less so compared to 
n = 2. The non-LTE line source function is not a Planck 
function, but instead is related to the ratio of the depar- 
ture coefficients (fe^'s) for the two levels involved in the 
transition. The ratio of the fe^'s for these three Balmer 
lines are shown in Fig. For all three Balmer lines the 
ratio of the &i's, and consequently the line source func- 
tion, decreases in the very region where the line cores 
form resulting in a fainter core compared to the wings. 
The LTE and non-LTE profiles are similar in the wings 
since the &i's approach unity at large optical depth. It is 
likely that other lines are also affected by non-LTE affects 
and the authors are working on a future paper that fully 
explores these issues for a number of atoms and ions. 

3.6. pre- CV Examples 
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Table 1. System properties for example pre-CVs 



Primary 



Secondary 



Object 


Teff(K) 


R(R©) 


M(M0) 


Toff(K) 


R(Rq) 


M{Mq) 


a(RQ) 


log (Fine) 


refs. 


GD 245 
NN Ser 
AA Dor 
UU Sge 


22,170 
55,000 
42,000 
85,000 


0.015 
0.019 
0.19 
0.34 


0.48 
0.57 
0.40 
0.63 


3560 
2900 
3000 
6250 


0.27 
0.17 
0.10 
0.54 


0.22 
0.12 
0.073 
0.29 


1.17 
0.95 
1.14 
2.46 


9.35 
11.32 
12.69 
13.79 


Schmidt et al. (199M 
jgatalan pt al. fl994) 

3ellet al. (l??!! 







The models and methods described above were appUed 
to four example systems (GD 245, NN Ser, AA Dor, and 
UU Sge) that have been well observed and documented in 
the literature. These systems were chosen because they 
represent the broadest range of properties for PCEBs, 
spanning 4 orders of magnitude in incident flux. Two 
systems, GD 245 and NN Ser, have WD primaries while 
the other two, AA Dor and UU Sge, have hot sub-dwarf 
stars as primaries. The purpose of this section is not 
to determine the properties of each system, but instead 
to illustrate some of the probable characteristics of their 
secondary's heated hemisphere. The orbital parameters 
and the bulk properties for the primary and secondary 
were taken from previous works and are listed in Table 

m 

3.6.1. GD 245 & NN Ser 

GD 245 and NN Ser are two examples of pre-CVs for 
which phase resolved spectroscopy has be en performed 
IjSchmidt et all 119951 iCatalan et allll994|) . Both sys- 
tems are observed to have Balmer lines composed of 
absorption wings from the primary and narrow emis- 
sion cores from the heated face of the secondary. The 
strength of the emission peaks at zero phase when the 
heated hemisphere is most visible. The observations of 
GD 245 clearly show the presence of additional emission 
lines from the secondary due to various metals (Mg II, 
Fe I, Na I, and Ca II). The spectra of NN Ser also show 
the presence of emission from Ca II and He I. Using the 
parameters listed in Table ^ niodcl atmospheres were 
computed for the secondary of both GD 245 and NN 
Ser. As described above, a single irradiated model was 
used to represent the average properties of the heated 
hemisphere. 

The temperature structures are shown in Fig. |H1 Both 
GD 245 and NN Ser have substantial temperature inver- 
sions at the top of the atmosphere and, in the case of 
NN Ser, extending down to ri.2 = 1. Both structures 
are compared to that of a non-irradiated M dwarf with 
\og{g) — 5.0 and Toff = 3500K (lowest dashed curve in 
Fig. 121). Apart from the temperature inversion, the ir- 
radiated model of GD 245 remains similar to a normal 
M dwarf atmosphere with convection present at layers 
where T1.2 < 0.1. However, the model for NN Ser, which 
included a much hotter primary (Tp = 55,000K) than 
GD 245, has a purely radiative atmosphere and is nearly 
isothermal down to T1.2 — 1. 

In the lower panel of Fig. |5| a total synthetic spec- 
trum is shown for the GD 245 system: i.e. the spec- 
trum of the WD primary and heated face of the sec- 
ondary were added using the parameters listed in Table 
n] assuming a 90° inclination. At A < 6OOOA, the pseudo- 
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Fig. 8 
irradiated models for the example cases listed in Table 
lines). Dashed lines are non-irradiated models with T;,ff 
(top) and Tcfi = 3000 A" (bottom). 



: 6000K 



continuum is shaped primarily by the WD's continuum 
and broad Balmer absorption lines. The red part of the 
spectrum is dominated by the flux from the secondary, 
which retains deep molecular absorption bands despite 
the extrinsic heating. In addition, the synthetic spec- 
trum approximately reproduced the basic collection of 
emission lines and relative line strengths shown in pre- 
viously published observed spectra (see Figs. 2 and 3 of 
Schmidt et al. 1995). There are certainly differences be- 
tween the observed spectrum and the synthetic spectrum 
(e.g., the number and widths of the metal lines). How- 
ever, given that only a single 1-D model was used for 
the secondary's heated hemisphere and only Hydrogen 
was treated in non-LTE, the s imilarities are encourag- 
ing. Also, Sch midt et ail l|1995l) estimated that GD 245 
emits about 60 times more H^ flux than predicted by 
simply counting the raw number of incident ionizing pho- 
tons from the WD. Schmidt et al.l l)1995D suggested that 
irradiation might reduce the number densities of molecu- 
lar and neutral atomic opacity sources and allow for the 
WD's Balmer continuum to supply the photons needed 
to explain the Hq flux. This idea is supported by the cur- 
rent models, which do show a dramatic decrease in the 
concentration of molecules and neutral atomic species. 
Furthermore, the synthetic spectrum shown in Fig. IHlhas 
Ha flux ~ 2 X 10"^'* ergs cm~^ sec~^ which is close to the 
observed value of 3.7 x 10~^^ ergs cm~^ sec~^ at phase 
zero. Note that a more realistic synthetic spectrum built 
from many irradiated models that simulate the variable 
heating across the secondary's surface reproduces the ob- 
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Fig. 9. — Combined spectra, i.e. primary flux plus the flux from the irradiated hemisphere of the secondary, for NN Ser (top) and GD 245 
(bottom). The atmospheric structures for the irradiated model atmospheres are shown in Fig. |S] The primary and s econdary spectra were 
co mbine d assuming an orbital inclination of 90° and phase 0. Note, both figures were designed to mimic those shown in ltjchmidt et all 119951) 
and lCat alan et al. (1994) to facilitate a comparison. Furthermore, both secondary spectra were rotationally broadened using Vrot = 70 km 
sec~^ and the resolution was reduced to 5 A for NN Ser an d 2A for GD 245 to matched the observed spectral resolutions from these studies. 
In the top panel, an observed spectrum of NN Ser from [Catalan efal. 11994) is also shown, but offset to make the plot more readable. 
Note, that the meanings of phases and 0.5 are reversed in the lCatalan et alj 119941) paper compared the present work. 



served Balmer line profiles of GD 245 very well l)Barmanl 
|2002}. 

The temperature structure of the heated hemisphere 
in NN Ser is significantly different from that of any non- 
irradiated main sequence star. The hot WD flux heated 
the secondary enough to suppress convection below the 
photosphere and reduced the presence of most molecules 
to very low concentrations. The combined WD + M 
dwarf synthetic spectrum for NN Ser is shown in Fig. 
El (top panel) assuming a 90° inclination. As with GD 
245, the flux from the hot primary determined most of 
the pseudo-continuum while narrow Balmer and Paschen 
lines filled in the cores of the WD hydrogen absorption 
lines. Also shown in Fig. |5|is a digitized version of the 
observ ed spectrum of NN Ser published bv lCatalan et al.l 
l) 1994ft for a similar phase. As with GD 245, the synthetic 
spectrum is very comparable to the observed spectrum of 
NN Ser. The pseudo-continuum and most emission lines 
are well reproduced by the model. The most noticeable 
differences are the strengths of the H emission lines which 
are weaker in the red half of the observed data. It should 
be stressed that no attempt has been made yet to find 
the best-fitting model for the observed spectrum. This 



comparison is simply for illustrative purposes and a more 
detailed analysis of observed spectra will be the subject 
of a later paper. 

3.6.2. AA Dor & UU Sge 

AA Dor and UU Sge are examples of extreme irradi- 
ation. The primaries, in both cases, are hot sub-dwarfs 
with radii more than an order of magnitude larger than 
the WD primaries of GD 245 and NN Ser. Furthermore, 
the primary of UU Sge is extremely hot, Tp = 85,000K. 
The atmospheric structures for the heated hemispheres 
of UU Sge and AA Dor are shown in Fig. |S1 Both mod- 
els have large temperature inversions and are radiative 
to depths well below the photosphere. Temperatures in 
AA Dor rise by nearly 20,000K between ri.2 = 1 and 
Ti.2 — 10~^. In the heated hemisphere model for UU 
Sge, temperatures are well above 20,000K with a maxi- 
mum temperature near ri.2 = 10~^ of almost 40,000K. 
The impact of this heating is also evident in their spec- 
tra. FigurelTUIshows the flux from the heated hemisphere 
model only; for these systems the dominant flux is from 
the primary and the secondary simply makes a small con- 
tribution to the overall pseudo-continuum. As in Fig. 
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Fig. 10. — Spectra for the irradiated hemispheres of A A Dor 
and UU Sge. The atmospheric structures for the irradiated model 
atmospheres are shown in Fig. |S] For both spectra, the resolution 
was reduced to 5A. 



a number of emission lines form in the bluest parts of the 
spectrum. 

Based on light curve analyses, it has been estimated 
that the temperature varies from 20,000K and 2000K 
across the heated a nd non-heated faces of AA Dor 
IjHilditch et alJ l2003j) . Such a difference is consistent 
with the irradiated and non-irradiated structures from 
Fig. IHl (keeping in mind that the models are averages 
for each hemisphere). The light curve analyses als o sug- 
gest that the albedo is very close to one (Hilditch et al.l 
l200.1ft which is also consistent with the model being fully 
radiative (and Fig. EJ. 

4. DISCUSSION AND CONCLUSIONS 

Model atmospheres have been computed for irradiated 
stars in pre-CV systems and high-resolution synthetic 
spectra are now available for direct comparison to ob- 
servations. The models presented above illustrate the 
dramatic differences that may exist between the heated 
and non-heated hemispheres of the secondary in pre-CV 
systems. Large temperature inversions, which can ex- 
tend down to the continuum forming layers of the pho- 
tosphere, are predicted on the heated hemisphere. Such 
temperature inversions naturally lead to emergent spec- 
tra that are unlike normal main sequence stars and are 
very different from black bodies. A steep Balmer-edge 
and strong hydrogen emission lines are present in addi- 
tion to a forest of narrow metal lines in emission. It has 
also been shown that departures from LTE are to be ex- 
pected and can have significant impact on the hydrogen 
line profiles. The present models are consistent with pre- 
vious findings for the atmospheric structures and bolo- 
metric reflection albedos (Vaz & Nordlund 1985; Nord- 
lund & Vaz 1990; BS93). 

The large differences between the irradiated and non- 
irradiated atmosphere models suggest that steep horizon- 
tal temperature and pressure gradients may exist across 



the surface of a pre-CV secondary. If present, the sec- 
ondary's heated hemisphere will have a SED that varies 
substantially across the heated surface, which may be 
underrepresented by single 1-D averaged models. Such 
gradients could also lead to strong horizontal flows and 
redistribution of en ergy to the non-irradiated hemisphere 
lIKirbivik fc SmithI 11973 IBeer fc Podsiadlowskil |2002). 
Furthermore, these gradients may lead to convectively 
unstable layers above the radiative-convective boundary 
predicted by the mixing length theory. The present work 
will be expanded to include full surface modeling of pre- 
CV secondaries in order to compute even more realistic 
light curves and synthetic spectra. 

A detailed non-LTE analysis is currently underway for 
a large number of atoms and ions for the strongly irra- 
diated systems. For the cooler cases that still have sig- 
nificant concentrations of molecules in their atmosphere, 
solving the rate equations is a very challenging problem. 
Much of the atomic and molecular data (e.g. coUisional 
cross-sections for atom/ion plus molecule interactions) 
are not yet available and are difficult to obtain. How- 
ever, some useful limits can be placed on the magnitude 
of the non-LTE effects for such cases. 

High-resolution, phase-resolved, spectra have been 
recently obtained for several interesting pre-CVs 
(Exter et al. 2003a y) and a future paper will be devoted 
to an analysis of these systems using synthetic spectra 
from irradiated model atmospheres. The example pre- 
CV systems described above are also prime candidates 
for detailed comparisons between observed spectra and 
synthetic spectra. In addition to providing insight into 
the atmospheric conditions of irradiated M dwarfs, such 
comparisons will also provide valuable tests for the mod- 
eling techniques that are currently being used in other 
areas (e.g. White Dwarf - Brown Dwarf pairs and close- 
in extrasolar planets). 
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